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THE SHAPE OF GEOLOGICA L HORTZOTsfg 

The invention relates to a process for locally 
5 determining the shape of geological horizons from data 
obtained by processing seismic measurements. 

The invention also relates to a device to determine 
- locally the shape of geological horizons from data obtained 
by processing seismic measurements. 
10 Finally, the invention relates to a computer software 

package permitting the operation of a programmable device 
according to the invention, for the practice of a process 
according to the invention. 

U.S. 6,151,555 relates to a process for processing 
15 seismic signals and a device for the determination of a cube 
of variance value. This process provides for computing for 
each node of three-dimensional seismic matrix a variance 
value equal to the variance of the seismic amplitudes 
adjacent the node in question. The vicinity of the node in 
20 question is a subset of the nodes of the three-dimensional 
matrix whose distance to the node in question is less than a 
selected threshold. This process relies on the hypothesis 
that the variance value thus defined is important for nodes 
near a discontinuity. 
25 The document WO 99/64896 discloses a process for 

interpretation of seismic data, particularly adapted to 
identify the limits of hydrocarbon reservoirs. 

A first object of the invention is to improve the 
analyses of the profile of the geometric horizons. 
30 A second object of the invention is to facilitate the 

determination of the unitary vectors normal to the horizons 
of a three-dimensional seismic cube. 
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A third object of the invention is to facilitate the 
determination of an index of curvature of the horizons in a 
three-dimensional seismic cube. 

The invention has for its object a process for locally 
5 determining the shape of geological horizons, in which there 
is constructured a continuous function Sij,k(t) by 
interpolation or approximation of the discrete seismic traces 
of a three-dimensional seismic matrix, said function being 
designated as "continuous local seismic trace", comprising 
10 the following steps: 

a) . using as optimum offset of two adjacent continuous 
local seismic traces, the offset value rendering maximal 
their correlation function; 

b) . taking as conditional neighborhood a "central" 
15 continuous local seismic trace Sij,k(t) , said neighborhood 

consisting in adjacent traces Spq,k(t) corresponding to 
optimum offsets associated with correlations Rij,pq,k(h) 
greater than a predetermined threshold comprised between 0 
and 1 ; 

20 c) . defining for each continuous local seismic trace 

Spq,k(t) of the conditional neighborhood a residual value 
relative to said "central" continuous local seismic trace 
Sij,k(t) comprising parametric coefficients; 

d) . determining the parametric coefficients by 

25 minimizing a set of residuals on the conditional 
neighborhood . 

According to other characteristics of the invention: 
- the residual values of step c) are absolute values of 
differential parametric shape, for example >c»y4>q,k«l%k.(p-0 + bij^«(q-j) - 
3 0 in which aij,k and bij,k are parameters and hij,pq,k is the 

optimum offset maintained between the "central" continuous 
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local seismic trace Sij,k(t) and an adjacent continuous local 
seismic trace Spq,k(t), 

minimizing the set of residuals Pij,pq,k on the 
conditional neighborhood by a minimization of a sum of powers 
5 of the residuals of the following type 

Prf in which a IS a power greater 

than or equal to 1, 



10 As the case may be, each term (pij,pq,k)*^ of this sum can 

be weighted by an increasing functional factor of the 
correlation coefficient Pij,pq,k (hij,pq,k) • 

- the coefficients aij^k and bij^k determined in step d) 
are used to define a unitary vector N(i,j,k) of coordinates 

15 N^'d^j^k), N^(i,j,k), N^(i,j,k), for example of the form: 



.... p'^" 



^ ° ^ mm = -j — ^ 

V(aiMc)'+(baj,)' + l 



25 and there is chosen this unitary vector as a "normal" 

unitary vector orthogonal to the horizon at the node (i, j, 
k) , 

- the set of normal unitary vectors N(i, j, k) is used 
to define the field of normal unitary vectors, for the 
30 display on a screen of the profiles of the seismic horizons, 
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there is selected as the index of curvature an 
increasing function of the residuals comprising the 
parametric coefficients determined in step d) . 

The invention also has for its object a device for 
5 practicing the process according to the invention, comprising 
means to use as optimum offset of two adjacent continuous 
local seismic traces, the offset value rendering maximal 
their correlation function, means to take as the conditional 
neighborhood of a "central" continuous local seismic trace 

10 Sij,k(t), the sub -neighborhood consisting in adjacent traces 
Spq,k(t) corresponding to optimum offsets associated with 
correlations Rij,pq,k(h) greater than a predetermined threshold 
comprised between 0 and 1, means to define for each 
continuous local seismic trace Spq,k(t) of the conditional 

15 neighborhood a value of residual relative to said "central" 
continuous local seismic trace Sij,k(t) comprising parametric 
coefficients, and means to determine the parametric 
coefficients by minimization of the set of residuals on the 
conditional neighborhood. 

2 0 According to other characteristics of the invention: 

- the device comprises memorization means and means for 
visualizing predetermined seismic parameters with the help of 
the process according to the invention. 

The invention finally has for its object a computer 
25 software package, comprising elements of program code to 
execute the steps of a process according to the invention, 
when said program is executed by a computer. 

The invention will be better understood from the 
description which follows, given by way of non-limiting 

3 0 example, with reference to the accompanying drawings, in 

which : 
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- Figure 1 represents schematically a three-dimensional 
seismic matrix and a continuous local seismic trace. 

- Figure 2 represents schematically an example of local 
neighborhood Nij^K consisting of a set of continuous local 

5 seismic traces Spq,k(t) adajcent the central reference trace 
Sij,k(t) itself located on the vertical of the geophone Gij and 
interpolating the seismic amplitudes about t=k. 

- Figure 3 represents schematically a vertical cross- 
section of conditional local neighborhood Nij,k(3r) of a 

10 continuous reference local seismic trace Sij,k(t:). 

- Figure 4 shows schematically a vertical cross- section 
of a three-dimensional seismic matrix with optimum offset 
hij,pq,k and a conditional neighborhood Nij,k(r) . 

- Figure 5 represents schematically a vertical cross- 
15 section of a three-dimensional seismic matrix analogous to 

Figure 4, with the presence of a fault or discontinuity. 

- Figure 6 represents schematically a functional diagram 
of a process according to the invention. 

With reference to Figures 1 to 5, identical or 
20 functional equivalent elements are designated by the same 
reference numerals . 

In Figure 1, a three-dimensional seismic matrix is 
obtained by reading measurements recorded by geophones Gij 
disposed on a network x,y at coordinate points i/j. The 
25 reading of the sampled measurements time-wise is represented 
along a descending axis t representative of the depth or of a 
vertical descent from the surface of the ground or the sea. 
The measurements are characterized by their amplitude, for 
example an amplitude read by the geophone G^^ at the time or 
3 0 at the sampling depth tk- The discrete measurement carried 
out by the geophone Gij as to time or depth tk is called the 
seismic amplitude Sij,k- 
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The set of seismic amplitudes corresponding to a 
geophone Gij of coordinates i,j is a uni -dimensional matrix 
(Siji, Sij2/ . . . . , Sijk/ . . . w SijN) called a discrete seismic 
trace, because this uni -dimensional matrix corresponds to the 
5 trace according to the point of horizontal coordinates i,j of 
the three-dimensional seismic matrix obtained by seismic 
measurements. 

The vertical axis t oriented along a descending 
vertical, usually designates time, but can also be considered 
10 as representing a depth from the surface. 

The invention also relates to the application of a third 
coordinate t representative of time, as well as a third 
coordinate t representative of depth. 

From the discrete seismic trace located on the vertical 
15 of a geophone Gij there is defined, by interpolation or 
approximation, discrete values about t=tk=k# a continuous 
function Sij,k(t) which is designated as "continuous local 
seismic trace". The methods of approximation or 

interpolation of discrete values to give rise to a continuous 
20 function are numerous, and comprise particularly polynomial 
interpolations or approximations, as well as polynomial 
trigonometric approximations or interpolations. 

Any other modification of interpolation or approximation 
providing a continuous function can also be used in the 
25 present invention to provide a "continuous local seismic 
trace" . 

In Figure 2, several continuous local seismic traces 
define an example of neighborhood of a "central" continuous 
local seismic trace Sij,k(t) as a reference. The neighborhood 
30 of a continuous local seismic trace Sij,k(t) is defined as the 
set of continuous local seismic traces whose horizontal 
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spatial indices p,q are adjacent horizontal spatial indices 
i,j of the reference continuous local seismic trace. 

By way of example, the horizontal spatial coordinates 
p,q corresponding to the geophone Gpq are adjacent horizontal 
spatial coordinates i, j corresponding to the geophone Gij if 
the absolute values of the differences i-p and j-q are less 
than given whole numbers, for example 2. 

In this case, as shown in Figure 2, the continuous local 
seismic trace Sij,ic(t) is associated with eight adjacent 
continuous local seismic traces surrounding the "central" 
continuous local seismic trace Sij,k(t). 

In the case of continuous local seismic traces produced 
by seismic measurements, the shape of the geological horizons 
introduces vertical offsets between adjacent continuous local 
seismic traces. So as to determine the relationships between 
two adjacent local continuous seismic traces Sij,k(t) and 
Spq,k(t) centered on the same vertical sampling coordinate 
t=tk=k and corresponding to different spatial coordinates i,j 
and p,q, there can be calculated the correlation function 
Rij^pq^k(h) of the two continuous local seismic traces. 

The correlation function of two adjacent continuous 
local seismic traces is obtained by the following formula 



in which the numerator corresponds to the covariance 
function of Sij,k(t) and Spq,k(t) obtained by the following 




expression 
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In this integral defining Cij,pq,k(h), the parameter A defines 
a "vertical investigation window" about t=tk=k. 

For example, if the continuous local seismic traces 
Sij,k(t) and Spq,k(t) are trigonometric polynomials of the 
5 following form, interpolating seismic data 

m 

10 

then it can be shown that the covariance function Cij,pq,k (h) 
is itself a trigonometric polynomial of the following form in 
which the coefficients Ae^^'^'^''' and Bs^^'^^''' depend on the 
coefficients as^^'"", as^"'''', ba^^''' and bs^*'''': 

A conventional and known mathematical result is that the 

2 0 correlation function R-gj»q.kW . / 

translates a similarity of correlated functions when this 
correlation function is near 1. 

The study of the correlation functions of the continuous 

25 local seismic traces Sij,k(t) and Spq,k(t) permits defining as 
optimum offset the value hij,pq,k/ of h corresponding to the 
maximum of the correlation function Rij,pq,k(h), which is to 
say to the maximum of the correlation function nearest 1. 
Contrary to the existing art based on a discrete formulation 

30 of the correlation function Rij,pq,k(h), the use of a 
continuous formulation of the seismic traces and hence of the 
correlation function permits obtaining optimum offset which 
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is not necessarily a whole multiple of the sampling step 
along the vertical axis corresponding to the variable t. 

The definition of the optimum offset hij,pq,k of two 
adjacent continuous local seismic traces permits obtaining a 
5 first approximation of the horizons passing through this 
reference continuous local seismic trace Sij,k(t), as shown in 
Figures 3 to 5 . Contrary to the prior art, the fact of using 
optimum offsets which are not required to be a whole multiple 
of the sampling interval along the vertical axis, permits 

10 avoiding errors known as "aliasing" . 

In these Figures 3 to 5, there is defined a conditional 
neighborhood Nij,k(r) of the reference local continuous 
seismic trace Sij,k(t) which is a sub -neighborhood of the 
initial neighborhood of the continuous local seismic traces 

15 used in the computation of correlation and of definition of 
the optimum offsets. The conditional neighborhood N^^'^(r) is 
selected such that, for any local continuous seismic trace 
Spq,k(t) belonging to Nij,k(r), the corresponding optimum offset 
hij,pq,k is such that the correlation Rij,pq,k (hij ,pq,k) between 

2 0 Sij,k(t) and Spq,k(t) is greater than a predetermined threshold 

r comprised between 0 and 1 . 

■ As is seen in Figure 5, the fixing of the threshold r 
also permits including in a conditional neighborhood a 
geological discontinuity or halt which constitutes an 
25 important advantage relative to the prior art. 

The mentioned arrangements thus provide a continuous 
modelization permitting use of the invention. 

In Figure 6, the digital modelization described in 
reference to Figures 1 to 5 is carried out in step 100 to use 

3 0 a process according to the invention. 

After this step 100 of digital modeling, there is 
defined in step 101, for each continuous local seismic trace 
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Spq,k(t) of the conditional neighborhood Nij,k(3r), a value of 
residual relative to said reference central continuous local 
seismic trace Sij,k(t) comprising parametric coefficients. 

The value of the selected residual Pij,pq,k translates the 
5 variation between the horizon at the vertical of the geophone 
Gpq and the plane tangent to the horizon at the vertical of 
the geophone Gij , 

By way of non-limiting example, there can be defined as 
the residual value an absolute value of differential 
10 parametric shape, of the type: 

/3»,pqj=-l aWc-(P-0+^ak ' in which aij^k and bij,k 
are parameters and hij,pq,k is the optimum offset maintained 
between the "central" continuous local seismic trace Sij,k(t) 
and an adjacent continuous local seismic trace Spg,k (t) . 
15 After having defined the set of residuals on the 

conditional neighborhood Nij^kCr), there are determined the 
parameters aij^k and bij,k which correspond to a minimization of 
the set of residuals Pij,pq,k on the conditional neighborhood 
in step 102 . 

20 The criteria of minimization selected is arbitrary, and 

can for example comprise a minimization of a sum of power 
greater than 1 of the residuals, of the form 

PA , m which n IS a real 

25 positive power greater than 1. 

As the case may be, each term (pij,pq,k)*^ of this sum can 
be weighted by an increasing functional factor of the 
coefficient of correlation Pij,pq,k (hij,pq,k) . 

The minimization carried out in step 102 produces as a 
30 result in step 103 the best estimation of the coefficients 
aij^k and bij,k which are directing coefficients of the horizon 
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passing through the point of horizontal coordinates i,j and 
of the time coordinate or depth coordinate k. 

In step 104, according to a first modification of the 
invention, there are determined the unitary vectors normal to 
the horizon passing through the spatial coordinate i, j and 
time coordinate or depth coordinate k by formulae of known 
type : 



The normal unitary vector thus obtained or any vector 
which is parallel to it, permits a visualization which is 
more effective than the visualizations of the prior art, so 
as to give seismic parameters and profiles of geometric 
horizons with improved contrast relative to the prior art. 

For another modification of the invention in step 105, 
there is selected as the index of horizon curvature at the 
point of spatial coordinates i, j and of the time coordinate 
or depth coordinate k, an increasing function of the 
residuals comprising the coefficients aij^k and bij^k determined 
in step 103 . 

There can particularly be used as the indicator or index 
of curvature the minimization criterion 
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used to compute the parametric coefficients aij^k and bij^k 
determined in step 103. 

There can also be used without departing from the scope 
of the present invention, any other increasing function of 
5 the residuals on the conditional neighborhood Nij^kC^r) . 

Thus, in the case in which the horizon is perfectly flat 
in the neighborhood or point of spatial coordinates i,j and 
of time or depth coordinate k, all increasing functions of 
the residuals will be equal to zero and the local index of 

10 curvature will as a result also be zero. 

On the contrary, in the case in which the horizon at the 
point of spatial coordinates i,j and of time or depth 
coordinate k is not flat, the index of curvature calculated 
in step 105 of the process according to the invention 

15 provides an indication or measurement of the difference 
between the horizon in question and the locally osculating 
plane at the point of spatial coordinates i,j and of time or 
depth coordinate k. 

The invention also permits by the computation of the 

20 index of curvature and its graphical representation on a 
screen, to improve the contrast of the seismic measurements 
provided by geophones in the form of a three-dimensional 
seismic matrix and thereby to avoid additional seismic 
measurements in the region of geological discontinuities or 

25 faults. 

The preferred practice of a process according to the 
invention consists in using a device for the practice of the 
process containing a computer software package executing the 
steps of the process according to the invention described 
30 with reference to Figures 1 to 6 . 

To this end, the device according to the invention is 
provided with suitable memorizing means to memorize and 
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compute the successive elements used in the process according 
to the invention. 

Preferably, the practice of the process according to the 
invention takes place at a work station for seismic 
5 computations by loading a software package according to the 
invention comprising elements of program code to execute the 
steps of the process according to the invention. 

The software package can comprise program code elements 
relating to steps 100, 101, 102, 103 and 104 described with 
10 reference to Figure 6. 

Alternatively, the software package according to the 
invention can comprise program codes elements to execute 
steps 100, 101, 102, 103 and 105 described with reference to 
Figure 6 . 

15 According to other modifications, a computer software 

package according to the invention can also comprise program 
code elements to execute at least the set of process steps 
described with reference to Figure 6. 

The invention described with reference to several 

2 0 embodiments is in no way thus limited, but covers on the 
contrary any modification of form and all variations of 
practice, including or not the presentation of intermediate 
results on a visualization screen. 
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